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Data for Drell-Yan (DY) processes on nuclei are currently available from fixed target experiments 
up to the highest energy of -^s — 40 GeV. The bulk of the data cover the range of short coherence 
length, where the amplitudes of the DY reaction on different nucleons do not interfere. In this 
regime, DY processes provide direct information about broadening of the transverse momentum of 
the projectile parton experiencing initial-state multiple interactions. We revise a previous analysis 
of data from the E772 experiment and perform a new analysis of broadening including data from 
the E866 experiment at Fermilab. We conclude that the observed broadening is about twice as large 
as the one found previously. This helps to settle controversies that arose from a comparison of the 
original determination of broadening with data from other experiments and reactions. 
PACS: 13.85.Qk; 24.85.-hp 
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INTRODUCTION 

Understanding how partons propagate in a nucleus 
reveals precious information about the early stages of 
hadronization that is difficult to obtain by other means. 
This understanding also plays a role in using partons as 
a sensitive probe of matter produced in relativistic heavy 
ion collisions. One of the important observables is broad- 
ening of the transverse momentum of produced particles 
due to multiple interactions of a parton in the medium. 
This effect can be studied for the nuclear medium in 
Drell-Yan reactions P, [3, S| and deep-inelastic scatter- 
ing 0, Q on nuclear targets. The former is an especially 
clean tool, since the produced dilepton, which has no 
final state interactions, carries undisturbed information 
on the transverse momentum of the parton initiating the 
reaction. 

A parton propagating through a nuclear medium ex- 
periences Brownian motion in the transverse momentum 
plane Q , with the increase of the mean transverse mo- 
mentum squared being proportional to the length of the 
path times the medium density p{z), 
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where z is the longitudinal coordinate. The energy de- 
pendent dimensionless factor C characterizes multiple in- 
teractions of the parton with bound nucleons. In the 
light-cone dipole description of the broadening process 
Hi) C is related to the small separation behavior of 
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The cross section aqq{rT, s) introduced in [lo| varies with 
transverse qq separation fx and energy. This quantity is 
difficult to calculate, since it involves nonperturbative ef- 
fects, but it can be fitted to data for the photoabsorption 
cross section and the proton structure function i^2(Q^, x), 
which probe a wide variety of separations and energies 
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Il2l |. Since each of the scatterings is accompanied 
by a loss of energy, there is a close connection between 
the multiple scattering that leads to momentum broad- 
ening and initial state gluonic energy loss [Hill. It was 
found, however, that this induced energy loss effects the 
DY process much less than the vacuum gluon radiation 
initiated by the interaction on the nuclear surface [l5| . 

In the parton model, a single interaction of a quark 
propagating through the medium can be viewed as a 
quark-gluon correlation in a nucleus in its infinite mo- 
mentum frame. The corresponding twist-4 term [16] can- 
not be calculated, but is rather fitted to data. This phe- 
nomenology has not been successful so far, with the pa- 
rameter extracted from data varying by a factor of 20, 
depending on assumptions [l^ and also demonstrating a 
strong process dependence [18 1. 

The first observation of nuclear enhancement of large 
transverse momenta was made in hadron-nucleus colli- 
sions [31 . A sensitive tool to measure the transverse 
momentum broadening of quarks propagating through 
nuclei is heavy dilepton production in the Drell-Yan re- 
action. The leptons have no interaction in the nuclear 
medium and thus carry undisturbed information about 
the transverse momentum of the quark. Nuclear broad- 
ening of the mean transverse momentum squared of the 
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dilepton is defined as 
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Perturbative QCD predicts dauY /dp^ (x (for trans- 
versely polarized virtual photons) at large px- There- 
fore, the mean value {p^) is divergent, and the broad- 
ening of Eq. ([3|) might be ill defined (see discussion in 
[20|). Of course, the experiment imposes a cut oflF. Thus, 
the experimentally observed broadening is finite, but it 
may depend on the experimental acceptance. This diver- 
gence, however, sometimes cancels in the broadening for 
several reasons. First of all, perturbative QCD predicts 
no nuclear effects at large pr, i.e., duDvip-^) / dp^ 
AdaDYipp)/dp^ (see, e.g., [2l|). Secondly, nuclear shad- 
owing is expected to vanish at large pt, especially at an 
energy as low as 800 GeV, the energy of our current anal- 
ysis. In this case, the divergent high-pT part of the inte- 
grations cancel in ((S]), and one can measure the nuclear 
broadening effect within a restricted interval ofpx where 
the data exhibit nuclear effects. 

An analysis of broadening in the Drell-Yan reaction 
was presented in Q using data from the E772 experiment 
at Fermilab. The results of the analysis for A{p'^)dy are 
depicted in Fig. [T] as function of A. 
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FIG. 1: Broadening of transverse momentum in production 
of heavy dileptons (squares), J/^ (circles) and T (triangles) 
in pA collisions at 800 GeV. The data are from the E772, 
E789, and E866 experiments at Fermilab 0, [H [H, HI] . The 
curves are results of a fit with the parametrization A(p|,) = 
D [{A/2y/^ - l]. 

Following Eq. ([1]), one can parametrize the A- 
dependence as 
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A fit to the data points for the Drell-Yan reaction in 
Fig. [1] results in D = 0.027 GeV^, which corresponds to 
C = 2in (P). 



The same figure presents nuclear broadening in J/^* 
and T production at the same energy. Comparing with 
this data and results available at lower energies, one 
can make the following observations: (i) the magni- 
tude of broadening in reactions of heavy dilepton and 
quarkonium production differ by a factor of 5, while one 
would expect the Casimir factor 9/4. Indeed, heavy 
quarkonia should be produced via gluons, which interact 
in the medium more strongly than quarks, (ii) The 
errors for Drell-Yan data are small compared to J/^ 
production, whereas the latter is measured with much 
higher statistics. (iii) Comparison with broadening 
measured at lower energies presented in Fig. [2] suggests 
that broadening in heavy quarkonium production should 
rise with energy, whereas it is seen to be rather fiat, or 
even perhaps drop, at high energies in the Drell-Yan 
reaction. 
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FIG. 2: Energy dependence of the parameter D in fitted 
to data for dilepton and heavy quarkonium production in pA 
and nA collisions. 



The latter observation is difficult to explain, since the 
cross section controlling the multiple interactions in the 
medium rises with energy, (iv) The dipole approach, 
Eq. (21), based on a phenomenological color dipole cross 
section fitted to data, predicts the factor C [8] to be 
about twice as large as one suggested by the data 
presented in Fig. [TJ 

In this paper we revise the result of the previous anal- 
ysis of the E772 data for Drell-Yan nuclear ratios, 
combining them with the brand-new data [25] from the 
E866 experiment for Drell-Yan cross section in pD colli- 
sions. The improvement in the analysis arises from the 
fact that the E866 experiment provides more accurate pD 
cross section data at small pr than that available from 
E772. Additionally, it was assumed with no good justi- 
fication that the same simple parametrization of the pt- 
dependence is valid for both pD and pA collisions, which 
dictated a specific shape for the ratios and led to a sub- 
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stantial underestimate of the errors for pT-broadening. 
Our new analysis is less model-biased, but it results in 
somewhat larger errors. We also perform here an anal- 
ysis of new data from the E866 experiment, which leads 
to a broadening about twice as large as in 3], although 
with larger experimental errors. This result resolves the 
problems stated earlier. 



LIGHT-CONE DIPOLE DESCRIPTION OF 
NUCLEAR BROADENING 

When a projectile scatters from a nucleus, one can 
treat the target nucleons as independent with some de- 
gree of accuracy if viewed in the rest frame of the nucleus, 
since the bound nucleons are sufficiently well separated 
in that frame. However, when boosted to the infinite 
momentum frame the same reaction looks different since 
the nucleons may start communicating via their parton 
clouds at small Bjorken x. Indeed the Lorentz contrac- 
tion is proportional to 1/x, therefore it is much stronger 
at large, than at small x. As a result, the parton clouds 
are still separated for x > Xc, but overlap and interact at 
smaller x, x < Xc- The border line Xc can be estimated 
as. 
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where Ra is the nuclear radius. Note that this is only 
an order of magnitude estimate. There is a numerical 
factor in ([5]) that varies considerably for different species 
of partons and the way the parton distribution is probed 



The interaction between partons at x < Xc leads to a 
modification of their transverse momentum distribution: 
suppression of small momenta, enhancement at interme- 
diate, and no effect at large transverse momenta. This 
phenomenon is known as the saturation effect ^1], or 
color glass condensate [2^. One should keep in mind, 
however, that onset of this effect happens only eA, x < Xc- 

The partonic interpretation of hard reactions at small 
Bjorken x is known to be Lorentz non- invariant. Deep- 
inelastic scattering (DIS) looks like absorption of the 
virtual photon by a quark (antiquark) that belongs to 
the target in the infinite momentum frame of the latter. 
The reaction looks, however, quite different in the tar- 
get rest frame, as the interaction with a qq fluctuation 
of the high-energy virtual photon. Of course, all observ- 
ables, including the structure function F2{x,Q^), which 
is proportional to the total photon-target cross section, 
are Lorentz invariant. 

The Drell-Yan reaction is usually interpreted as an- 
nihilation of a quark with an antiquark belonging to the 
beam and target respectively (or vice versa) , qq II [30| . 
While this is correct in the rest frame of the dilepton, this 



interpretation of a heavy dilepton production is not ap- 
plicable in target rest frame. Indeed, in order to satisfy 
the kinematics annihilating with the projectile, the target 
parton must move with momentum ~ Pt/3^2j where 
X2 in standard notation is the Bjorken x of the target 
parton. At small X2, X2 <C 1, this would mean presence 
of highly energetic partons in the stationary target. The 
proper space-time interpretation in this case is analogous 
to the Weitzacker- Williams mechanism of electromag- 
netic bremsstrahlung. Namely, a beam quark (antiquark) 
develops a fluctuation, q qj* (with 7* ^11), which in- 
teracts with the target freeing the dileptonjsj, (s^]. 

What looked like overlap of parton clouds at small X2 
in the infinite momentum frame, now takes the form of 
interference between the amplitudes of a parton interact- 
ing with different nucleons. On the same footing it can 
also be interpreted as the lifetime of the \q'y*) fluctuation, 
given by the uncertainty principle 
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where gi, is the longitudinal momentum transfer in the 
q — > qj* fluctuation. The effective mass squared of the 
fluctuation is M^^. = M^,/a + my{l~a) + k'^/a{l-a), 
where a is the fraction of the light-cone momentum of the 
parent quark taken by the dilepton (7*), and fcy is the 
7* — 9 relative transverse momentum. Hereafter we do 
not discriminate between coherence time and length, as- 
suming that the quark has the speed of light. Evidently, 
the same length ~ l/^L defines the maximal longitu- 
dinal distance between two scattering centers that can 
interfere. This is why the distance, Eq. is also called 
coherence length. 

In the case of a nuclear target, the coherence length is 
of great importance. If this time interval is small com- 
pared to the mean spacing of nucleons in the nucleus, 
dileptons are radiated in the Bethe-Heitler regime, i.e., 
with no interference between the amplitudes of radia- 
tion from different nucleons. At the another extreme, 
Ic 3> Ra, the Landau-Pomeranchuk phenomenon governs 
the radiation. Namely, interferences suppress radiation 
with small transverse momenta, enhance it at intermedi- 
ate values of few GeV, and make no changes at higher 
momentum transfer. In both cases, the linear depen- 
dence on the path length, or nuclear thickness, Eq. ([T]), 
holds, although with different factors C 20]. The bor- 
derline between the two regimes, Ic = Ra, corresponds 
to the relation for Xc, Eq. ([5]). 



Short coherence length 

The mechanism of nuclear broadening depends on how 
long the coherence length is compared to the nuclear size. 
In the regime of short coherence length, Ic ^ Ra, the 
dilepton fluctuation appears for a very short time deep 
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inside the nucleus and is immediately released on mass 
shell by the interaction, which is assumed to be the same 
as on a free nucleon. Therefore, this interaction does 
not lead to any change in the transverse momentum de- 
pendence, i.e., it provides no broadening. However, the 
initial state interactions of the projectile quark affect the 
transverse momentum distribution of the quark, which 
acquires the shape W[ 
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Here r2|„ {fi , ^2 ) is the density matrix describing the im- 
pact parameter distribution of the quark in the incident 
hadron. We assume a Gaussian shape, 
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where (pq) is the mean value of the primordial transverse 
momentum squared of the quark. 

The transverse momentum of the radiated dilepton re- 
flects the broadening acquired by the quark in the nu- 
clear medium. However, the effect is weaker than for the 
quark, since the transverse momenta of the dilepton and 
quark are connected by the relation. 
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where a < 1 is the fractional momentum of the II. This 
relation clearly demonstrates that the effect of broad- 
ening cannot be translated into a modification of the 
quark distribution function (as is frequently done), in 
the regime of short coherence length. Thus, it confirms 
the breakdown of fcT-factorization imposed by the initial 
state interaction, as suggested in 33 1. 



No-shadowing sum rule 

To the extent that the coherence length for the Drell- 
Yan reaction is short compared to the internuclcon spac- 
ing in the nucleus, no shadowing is possible. This means 
that in spite of the nuclear modification of the transverse 
momentum dependence of Drell-Yan cross section, the 
number of quarks is conserved. This fact can be repre- 
sented in the form of a sum rule. 
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The nuclear ratios measured in the E772 and E866 ex- 
periments at SOOGeV depicted in Fig. [3] have a tendency 
to satisfy this constraint. Indeed, the ratios are below 
one at small pT, but rise above one at larger momen- 
tum transfers. The weak nuclear suppression observed 
at small pr, where the absolute cross section is maximal. 
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FIG. 3: The nucleus-to-deuterium and -beryllium ratios of 
dilepton production cross sections as function of transverse 
momentum. The data are from the E772 0] and E866 ^ 
experiments at 800 GeV. 



seems to be compensated by a stronger enhancement at 
larger pT in accordance with the sum rule Eq. (fTO|) . Data 
in the small px < 1 — 2 GeV domain, where the main 
part of the cross section comes from, are most important 
for broadening. In this region the E866 data in Fig. 3 
demonstrate stronger nuclear suppression than the E772 
data. The observed deviation from unity at small px, 
which is responsible for broadening, is about twice as 
large in E866 as in E772 data. 

Notice that the pT-dependence of the ratios in Fig. [3] 
has the typical shape predicted for the regime of satura- 
tion (of quarks), or color glass condensate [i^. There- 
fore, it might be tempting to interpret such data as an 
observation of this phenomenon. We warn again that 
this would be a wrong conclusion, since coherence and 
saturation are impossible in the short coherence length 
regime; initial state interaction effects imitate saturation. 



Long coherence length 

In this regime, Ic ^ Ra, the amplitudes 7* II for 
radiating a heavy photon from different nucleons interfere 
with each other. This interference, which modifies the 
Pt distribution of radiation, is known as the Landau- 
Pomeranchuk effect and is a genuine effect of saturation. 

Alternatively, one may say that a long-lived fluctua- 
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tion q —> j*q that propagates through the entire nucleus 
may be freed as a resuh of the multiple interactions, in 
contrast to what happens in the regime of short Ic, where 
the fluctuation is freed after an interaction with a single 
nucleon. These multiple interactions also lead to a dif- 
ferent type of broadening mechanism. Its origin can be 
understood from the fact that most of the q — > 7*5 fluc- 
tuations in the quark are invisible in the sense that their 
interactions with the target are so weak that the target 
has restricted resolving power or ability to discriminate a 
hard 57* fluctuation, with large intrinsic transverse mo- 
mentum kx , from a bare quark. However, because of the 
larger momentum transfer allowed by the multiple inter- 
actions in the nuclear medium, the resolution provided by 
the nucleus is greater than that of a free nucleon. The 
nucleus can therefore free harder fluctuations than a free 
nucleon. This is the source of nuclear broadening in the 
large coherence length regime. 

The same effect can be interpreted within the color 
dipole approach in terms of color flltering. The nucleus 
filters out projectile dipoles of smaller size than is pos- 
sible on a free nucleon, leading to larger transverse mo- 
menta. In this case, instead of Eq. ([7]) , the transverse mo- 
mentum distribution of dileptons radiated by a quark in- 
teractin g wi th a nucleus is given by the factorized dipole 
formula 13211 . 



da(qA IIX) 
d^qx 

X *t_^.(fi,a)4'„-(f2,a) 



iqT(ri-r2) 



1 - ^-h''^,(ri,X2)TA(b) 



Here ^q^- (r, a) is the pQCD calculated light-cone dis- 
tribution amplitude, which describes the dependence of 
the I (77*) Fock component on transverse separation and 
fractional momentum. The nuclear thickness function 
at impact parameter b is given by the integral of the 
nuclear density along the quark trajectory, TA{b) = 
PAih, z). 

The ratio of the nuclear to nucleon cross sections of the 
Drell-Yan reaction calculated in 32] has a shape typical 
for the Cronin effect. Such a modification of the trans- 
verse momentum dependence is a genuine effect of satu- 
ration (or its onset). In spite of the similarity with the 
analogous results in the short Ic regime, the mechanisms 
and their formulations in Eqs. ((7|)- (fTT|) are different. Ad- 
ditionally, the magnitude of the Cronin effect in the long 
Ic regime is smaller poj . 



ANALYSIS OF E772/866 DATA FOR DRELL-YAN 
REACTION 



Here we determine the nuclear broadening A(py) of 
dileptons produced on different nuclear targets relying 



on data for pj^-dependent ratios depicted in Fig. [31 
These include the nucleus-to-deuterium ratios R^^^ from 
the E772 0] experiment, and nucleus-to-beryllium ratios 
jlA/Be -fj.Qj^ ^j^g gggg experiment Q, both at 800 GeV. 

These results are quite reliable since they are least af- 
fected by systematic uncertainties. The ratios Ri are 
available in 7 bins in pT, i = I — 7. Correspondingly, one 
can calculate the broadening, Eq. ([3]), as 



T.^(^^RiPi Y^i^iPi 



where for each i-th bin 
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with pf the mean value of p^ within the i-th bin. In 
addition to the seven pr-intervals in data depicted in 
FiglSl we have included an eighth interval in Eq. ([72]) 
covering (pt)™"^ < pt < 00 and assumed that i?g = 1 
in accordance with perturbative QCD expectations (2lj 
(see, however, next section). Since the cross section falls 
steeply at large px, this point should not substantially 
affect the results. 

The differential cross section for the Drell-Yan reaction 
on a free nucleon, da^y/d'^PT, is more affected by sys- 
tematic uncertainties than the ratio is. We rely on the 
recently published results of the E866 experiment for the 
differential cross section in pp and pD collisions 25| . We 
use this data in for the analyses of all data presented 
in Fig. [3 

We parametrize the differential cross section of the 
Drell-Yan reaction as 



da 
dpip 



= N 



(14) 



The factor in the numerator is introduced to describe a 
possible forward minimum in the cross section indicated 
by data 

We found that the shape of the pr dependence does 
not vary with dilepton mass and X2^ as one can see in 
Fig. m A global fit to Drell-Yan data from pD colli- 
sions measured by the E866 experiment for different M 
and X2 bins [25| with common parameters n, and dif- 
ferent normalization factors N , led to quite a good de- 
scription, = 83.3 with Nd.o.f. = 52. The values of 
the parameters are Ai = 0.78 ± 0.08; A2 = 0.65 ± 0.04; 
A3 = 3.31 ± 0.21; n = 7.00 ± 0.66. The results of the 
fit and the data are depicted in Fig. H) The error bars 
include both statistical and systematic errors added in 
quadrature. 

Note that at large pr the pD cross section in Eq. (fT4[) 
is falling asp^^"*, much steeper than one may expect from 
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FIG. 4: Differential cross section of Drell-Yan reaction in 
proton-deuteron collisions. Data are from the E866 exper- 
iment [sJl at 800 GeV for different intervals of dilepton ef- 
fective mass. Curves show the result of the global fit with 
common parameters, Eq. (|14|l . controlling the shape of the 
cross section. 



FIG. 5: Results of the present analysis for broadening in the 
Drell-Yan reaction on different nuclei. Closed squares corre- 
spond to E772 data for C, Ca, Fe and W. Open squares 
correspond to E866 data for Fe and W. For comparison, 
broadening for J/^ and T, the same as in Fig. [T] are also 
shown. 



perturbative QCD. This happens due to restrictions im- 
posed by finiteness of energy, namely, both a;iand X2 rise 
with pt towards the kinematic limit causing an addi- 
tional suppression of the cross section [20]. Addition- 
ally, available data cover a rather restricted interval of 
Pt < 4 GeV. 

With the cross section in Eq. we performed an 
analysis of data from both experiments, E772 and E866. 
While the ratios measured in E772 are free from major 
systematic uncertainties, the absolute cross sections are 
substantially less reliable. Indeed, comparison of the dif- 
ferential cross sections of the Drell-Yan process measured 
for p ~ D collisions in the two experiments reveals sub- 
stantial differences. Therefore, we rely on the p — D cross 
section from the E866 experiment. The results of such 
a combined analysis of the E772 ratios are depicted in 
Fig. [5] by closed squares. 

The ratios measured in the E866 experiment are for a 
nucleus relative to beryllium. We use the chain relation 
between the two ratios, 



Ra/d — Ra/Bc Rbc/d 



(15) 



however, the beryllium-to-deuterium ratio is unfortu- 
nately lacking. Nevertheless, this cannot be a significant 
correction, since beryllium is a light nucleus. Indeed, the 



carbon to deuterium ratio depicted in Fig.[3]is compatible 
with no nuclear effects at all. If we simply fix Rbc/d — 1; 
the broadening given in Eq. (jl2p will be underestimated. 
However, one can do better by relying on a theoretically 
predicted value for Rse/D |26| with Eq. ([7]). Of course, 
this introduces an uncertainty related to model depen- 
dence, but this uncertainty should not exceed few percent 
of the observed broadening. The results of this analysis 
of the E866 data are depicted in Fig.[S]by empty squares. 

Notice that the cross section ratios are affected by sys- 
tematic uncertaities much less than the absolute cross 
sections. The ratios measured in both experiments have 
mainly an overall systematic uncertainty 1% in E866 3J], 
and 2% in E772 [2] measurements. One can see from 
Eq. (|12p that an overall variation of all Ri leaves the 
broadening unchanged. Therefore we could disregard the 
systematic errors of the ratios. 

If we fit the data in Fig.[5]assuming a linear dependence 
on A^/'^, in accordance with ([4]), and do this separately 
for the E772 and E866 data, we obtain. 



D{E772) = (0.029 ± 0.008) ( GeY/cf : 
D{E866) = (0.059 ± 0.009) ( GeY /cf 



(16) 



The value of D{E772) is close to the result of previous 
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analysis [Sj] depicted in Fig. [U although with a consid- 
erably larger error. However, the E866 data suggest a 
much stronger broadening. Nevertheless, one may con- 
sider the results Eq. (fTB|) as nearly consistent since they 
lie within two standard deviations. If the broadening ac- 
cording to the E866 data were actually about twice as 
large as that given by the previous analysis 3] of the 
E772 data, this would settle the contradiction with theo- 
retical expectations and conclusions drawn from the data 
for the broadening of J/^* and T mentioned above. 



BEYOND THE DATA 

Although QCD predicts no nuclear effects at high px, 
the data in Fig. [3] show considerable nuclear enhance- 
ment at intermediate values of px- This is the well- 
known Cronin effect first observed for hadronic beams 
[l9l|. It is not clear from the data how far in px this 
effect can propagate. The assumption made above for 
the last unobserved pr interval, i?8 = 1, is probably in- 
correct, therefore the broadening effect shown in Fig. [5] 
could be be underestimated. The ad hoc shape of R{pt) 
used in Q is difficult to justify and docs not suggest a 
solution to the problem. 

In order to get a hint about how much of the effect 
might have been missed, one may make use of a theo- 
retically calculated value for R^. We rely on the dipole 
formalism, Eqs. ([T])-®, from which one is able to cal- 
culate the pt dependence of the Drell-Yan reaction on 
proton and nuclear targets in a parameter-free way and 
explain the data in Fig. [3] well with no adjustment of 
parameters. The dipole approach also provided the first 
quantitative explanation of the Cronin effect in available 
data for hadron-nucleus collisions and correctly predicted 
the effect for RHIC 21]. Therefore we believe that a cal- 
culated value of i?8 should be rather reliable. Doing this, 
Eq. dm leads to somewhat larger values for the broad- 
ening. 

The fit of Eq. (jj]) to this data leads to new values of 
parameter £), 



BROADENING VERSUS COHERENCE LENGTH 

The mechanism of broadening is expected to correlate 
with the coherence length [1, llo, 21 1 as described above. 
The coherence length of the Drell-Yan process is con- 
trolled by the value of X2, Ic ^ 1/ 171^x2 [15]. Thus, at 
large X2 the dilepton is produced inside the nucleus and 
only half of the nuclear thickness contributes to broaden- 
ing. However, at small X2 the Drell-Yan pair is radiated 
outside the nucleus and the entire nuclear thickness is 
effective. It would be a spectacular observation to see a 
manifestation of this effect in data. 

Using recent data from the E866 experiment for the 

we can also test 



PT-dependent ratios binned in X2 |34l . 
the predicted a:2-independence of broadening. We have 
performed a similar analysis of the Fe/Be and W/Be 
ratios in different X2 bins and the results are depicted in 
Fig.lH 



> 
O 



0.5 



0.4 



0.3 



0.2 



0.1 



: Fe/Be 
- * W/Be 












c 


C 


c 


) 


c 










00 


0.05 




0.1 



X2 

FIG. 6: Broadening of the mean transverse momentum 
squared in the Drell-Yan reaction measured [s^, [s^ on tung- 
sten and iron targets relative to beryllium, A(py) = {pj')a — 

{PT}Be. 

Unfortunately, within the rather large statistical errors 
one cannot draw a definite conclusion regarding a possi- 
ble variation of broadening with X2 ■ 



D{E772) = (0.038 ± 0.008) ( GeY/cf ; 
DiE866) = (0.069 ± 0.009) ( GeY/cf . (17) 

This procedure involving theoretical calculations may 
be considered as a means to estimate the systematic un- 
certainty. The sign (positive) we obtain for the uncer- 
tainty is certainly correct, and its magnitude is not large. 
Note also that these corrections are correlated for the 
two sets of data, and therefore one must apply them con- 
sistently to both the the results of the E772 and E866 
experiments. 



CONCLUSIONS 

Analysis of data from the E866 experiment for the DY 
reaction on nuclear targets has resulted in a determina- 
tion of the magnitude of the nuclear broadening parame- 
ter _D(£'866) given in ^TE\i and ^T7\ . which is about twice 
as large as the one suggested by the E772 data. This 
result helps to resolve the controversies that were a par- 
tial motivation for this analysis. One of these was the 
expectation that the observed broadening in DY would 
be a factor of 4/3 less than that observed in J/\l/ and T, 
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as dictated by the Casimir factor. The updated result is 
consistent with this expectation and also with the expec- 
tation of a rising energy dependence of broadening. And, 
last but not least, the value of D{E866) agrees with the 
parameter-free theoretical prediction of ^] . 

This new result also agrees with broadening observed 
by the HERMES experiment at HERA in hadron pro- 
duction in deep-inelastic scattering on nuclei fij . Indeed, 
both measurements are well described within the same 
color-dipole approach [s^. 

Since most of the data from both experiments E772 
and E866 correspond to the regime of short coherence 
length, the observed broadening is directly connected via 
the relation in Eq. ([9]) to the broadening of a quark propa- 
gating through the nuclear medium. In the regime of long 
coherence length the mechanism of broadening changes 
and takes the form of color filtering in hard coherent 
scatterings. To see a possible variation of broadening as 
function of coherence length, we performed an analysis 
of data additionally binned in X2- Unfortunately, the 
restricted statistics did not allow us to draw a definite 
conclusion from this analysis. 

Broadening of the transverse momentum of a parton 
in a medium leads to induced gluon radiation, i.e., to 
induced energy loss. This is the basis of the effect called 
jet quenching which is considered to be a sensitive probe 
for properties of matter created in heavy ion collisions. 

Note that broadening itself is much less dependent on 
models for hadronization, than jet quenching. Therefore 
it is a better and more certain probe. 

Acknowledgments: This work was supported in part 
by U.S. Department of Energy, the Research Ring "Cen- 
ter of Subatomic Studies" (Chile), Fondecyt (Chile) 
grants 1030355, 1050519 and 1050589, and by DEC (Ger- 
many) grant PI182/3-1. 



[1] 
[2] 
[3] 

[4] 
[5] 

[6] 
[7] 
[8] 



NAIO Collaboration, P. Bordalo et al., Phys. Lett. B193, 
373 (1987). 

E772 Collaboration, E772 Collaboration, D.M. Aide et 

al., Phys. Rev. Lett. 66, 2285 (1991). 

P.L. McGaughey, J.M. Moss, and J.-Ch. Peng, NU- 

COLEX 99, Wako, Japan, 1999 (h ep-ph/9905447 ); Ann. 

Rev. Nucl. Part. Sci. 49, 217 (19997 

HERMES Collaboration, A. Airapetian et al., Eur. Phys. 

J. C20, 479 (2001). 

W. Brooks, talk presented at the Workshop "Par- 
ton Propagation through Strongly Interacting Systems" , 
ECT*, Trento, 26 September - 7 October, 2005. 
J. Hufner, Y. Kurihara, H.J. Pirner, Phys. Lett. B215, 
218 (1988). 

J. Dolejsi, J. Hiifner and B.Z. Kopeliovich, Phys. Lett. 
B312, 235 (1993). 

M.B. Johnson, B.Z. Kopeliovich and A.V. Tarasov, Phys. 
Rev. C 63, 035203 (2001). 



[9: 

[10 

[11 

[12; 

[13 
[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 

[21 

[22: 

[23: 

[24: 

[25: 

[26: 

[27 

[28: 

[29 
[30 
[31 



[32: 



[33 

[34' 
[35: 

[36] 



A. Accardi and D. Treleani, Phys. Rev. D64, 116004 

(2001) . 

B. Z. Kopeliovich, L.I. Lapidus, and A.B. Zamolodchikov, 
Sov. Phys. JETP Lett. 33, 595 (1981); Pisma v Zh. Ex- 
per. Teor. Fiz. 33, 612 (1981). 

K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D59 
014017(1999). 

B.Z. Kopeliovich, A. Schafer and A.V. Tarasov, Phys. 
Rev. D62 054022 (2000). 

R. Baier, Yu.L. Dokshitzer, A.H. Mueller, S. Peigne and 

D. SchiflF, Nucl. Phys. B 484, 265 (1997). 

I. Vitev, Phys. Lett. B562, (2003) 36. 

M.B. Johnson et al. [FNAL E772 Collaboration], Phys. 

Rev. Lett. 86, 4483 (2001); Phys.Rev. C65, 025203 

(2002) . 

M. Luo, J.W. Qiu, and G. Sterman, Phys. Rev. D50, 
1951 (1994). 

J. Qiu, X. Zhang, Phys. Rev. Lett. 84 5049 



X. Guo, 
(2000). 

E. Wang and X.-N. Wang, Phys. Rev. Lett. 89 162301 
(2002). 

D. Antreasyan et al., Phys. Rev. D 19, 764 (1979). 

B. Z. Kopeliovich, J. Raufeisen, A. V. Tarasov and 

M. B. Johnson, Phys.Rev. C67, 014903 (2003). 

B.Z. Kopeliovich, J. Nemchik, A. Schafer and 

A. V. Tarasov, Phys. Rev. Lett. 88, 232303 (2002). 
E772 Collaboration, D.M. Aide et al., Phys. Rev. Lett. 
64, 2479 (1990); M. Leitch et all., Nucl. Phys. A522, 
351c (1991). 

E789 Collaboration, M.J. Leitch et al., Phys. Rev. D52, 
4251 (1995). 

E866 Collaboration, M.J. Leitch et al., Phys. Rev. Lett. 
84, 3256 (2000). 

E866 Collaboration, J.C. Webb et al., |hep-ex/0302019[ 
to appear in Phys. Rev. Lett. 

M.B. Johnson, B.Z. Kopeliovich, and I. Schmidt, 
'hep-ph/0701015 

B. Z. Kopeliovich, J. Raufeisen, and A.V. Tarasov, Phys. 
Rev. C62, 035204 (2000). 

L.V. Gribov, E.M. Levin and M.G. Ryskin, Nucl. Phys. 
B 188, 555 (1981); Phys. Rep. 100,1 (1983). 
L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 
(1994); D 49, 3352 (1994). 

S. D. Drell and T. Van, Phys. Rev. Lett. 25, 316 (1970) 
[Erratum-ibid. 25, 902 (1970)]. 

B. Z. Kopeliovich, proc. of the workshop Hirschegg 
'95: Dynamical Properties of Hadrons in Nuclear Mat- 
ter, Hirschegg January 16-21, 1995, ed. by H. Feld- 
meyer and W. Norenberg, Darmstadt, 1995, p. 102 
(|hep-ph/9609385). 

B.Z. Kopeliovich, A. Schafer and A.V. Tarasov, Phys. 
Rev. C 59, 1609 (1999) (see extended version in 
hep-ph/9808378). 

G.T. Bodwin, S.J. Brodsky, and G.P. Lepage, Phys. Rev. 
D 39, 3287 (1989). 

M.A. Vasiliev et al., Phys. Rev. Lett. 83, 2304 (1999). 
See plots and tables in the E866 web page, 
http: / /p25ext .lanl.gov /e866 /papers /papers. html 
B.Z. Kopeliovich, J. Nemchik, E. Predazzi, 
A. Hayashigaki, Nucl. Phys. A740, 211 (2004). 



